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The wall interference is obtained for a wind tunnel of 
elliptic section for the tws* cases of closed and open work- 
ing sections. The approximate and exact methods used gave 
results in practically good agreement. Corresponding to 
the result given by Glauert for the case of the closed rec- 
tangular section, the interference is found to be a minimum 
for a ratio of minor to major axis of 1: J2. This, however, 
is true only for the case where the span of the airfoil is 
small in comparison with the width of the tunnel. For a 
longer airfoil the favorable ellipse is flatter. In the 
case of the open working section the circular shape gives 
the minimum interference. 



INTRODUCTION 



The wall interference exerted on a model during wind 
tunnol tests has been given, for the case of a circular 
boundary, by Prandtl (reference l), for the closed rectan- 
gular boundary approximately by Glauert (reference 2), and 
accurately by Terazawa (reference 3). Recently Theodorsen 
(reference 4) combining ^the cases of open and closed sec- 
tions computed the wall— int erf er ence coefficient for the 
rectangular wind tunnel. However, apart from the tests of 
Knight and Harris (reference 5) on the expected wall in- 
terference of the projected large— scale wind tunnel of 
elliptic section, there have as yet appeared no investiga- 
tions on the elliptic section. Whether for a large— scale 
tunnel or one of usual size the determination of the wind- 
tunnel correction for free-air conditions is an important 
problem. 

•Journal of the Society of Mechanical Engineers, Japan, 
vol, 36, no, 190, Feb. .1933, pp. 123-127. 
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The problem of the interference exerted by the tunnel 
walls is. treated with the aid of the known results of the 
vortex— filament thoery applied to the boundary. In the 
present paper the wall interference is treated as a Dirichlet 
problem. Depending on whether the "boundary conditions are 
completely satisfied or satisfied at only a determined num- 
ber of points an exact or approximate solution, respectively, 
is obtained. The reason. that both are given is to show that 
for practical purposes the difference in the results obtained 
is negligible and to emphasize the usefulness of the latter 
method when the former is not available. 

Ihe closed and open working sections are characterized 
by the conditions that the normal and tangential components 
of the velocity at the wall are, respectively, zero. The 
former is represented by the N.P.L. wind tunnel, the latter 
by the G-ottingen and Eiffel tunnels. 

APPROXIMATE METHOD POR THE CLOSED WORKING-' SECTION 



Pigure 1 shows the two symmetrical vortices for the 
case of the closed section of elliptic boundary. The com- 
plex potential f = cp + ity of the flow within the ellipse 
is divided into two parts: f x due to the effect of the 
vortex alone, and . f s due to the existence of the wall. 
Por fj_ it is clear that 

f., is as yet unknown, but since it must be an analytic com- 
plex function in the field under cons iderat ion, put 

as 



+ / (a n + *Pn> * n < 2 > 



. f 3 = a Q + 19 

i i 

n=l 

where a n , 0 n ( n = 0 » 1»2)> are real con-stants. Introducing 
polar coordinates (r, 8) with the aid of the relation 

z = ro^ tho velocity components in the directions of r 
and 0 from (l) and (2) are obtained, respectively, as 
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v - * - 



£s (r 2 +s 2 )sin 9 

n ( r 2 -8 2 ) 2 +4 r 2 s 2 b in 2 9 



CO 



n = l 



- IM - Us. (r s -s a )cos9 
9 " r be ~ rr ( r 2 -s 2 ) 8 +4r 2 s 8 sin 2 9 



>► (3) 



CO 



'I 



nr n "" l (ons in n9 + 3 n cosn6) 

Henca, if X is the angle made "by the normal of the ellipse 
with the radius vecter, the normal component of the velocity 
at the boundary is 



r H " [vr] r=CTa cosX + [v e ] r=aa 



s in X 



where 



Ta S a sin ( 9-X) + k 2 a in (9 + X ) 
rra 8 ' (ff 2 ~k s ) a + 4cr s k 2 s in s 9 

co 

+ y na n_l < ,. il a n - l co3(n9 + X) 

CO 

• j^T na n "" l p n o n '" X *in(n 8 + x ) 

n = l 



(4) 



k = a/a, a = [ ( 1- 1- e 8 cos 2 9 ) ] 

tan X = € (l-c^) l/a sin 29/2(l~e 2 cos s 9) 
e = (l-D 2 /a s ) l/a = eccentricity 



(5) 



since there is no normal velocity component at the "boundary 
for the closed working section vjj = 0, so that from this 
condition the constants a n , & n are determined. Prom the 
condition" of symmetry there is immediately ot n = 0 and, more- 
over, p n vanishes for even n, Hence writing 
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na n;_1 3n=- rsB n /na 3 (6) 
gives the equation for determining the constants 

B 1 sin(e+X)+ B 3 CT S sin(38 + X) + B 5 cr*sin(58 + \) + ... 

= 0 g sin (8 - \) + k a sin(8 + \) . ? ^ 

C<j s - k 8 ) 2 + 4a s k 2 sin :3 8 

In this relation 8 takes all values from 0 ... , to 2tt. Since 
\ and a are functions of 8 , .it is difficult to solve 
this relation accurately. However, for special values of 8 
satisfying this condition a relatively easy approximate com- 
putation can be made. Since, on account of the symmetry, 
it is sufficient to consider only one quadrant of figure 1, 
it is possible to choose, for example, for 8 the four 
values 15°, 30°, 60°, and 90°, keep four terms of series on 
the left side of (?) and obtain four independent first— degree 
equations from which to solve for B lf B a , B 5 , and B 7 . Since 

f s is not a very rapidly convergent series the number of 
equations may be large. It should be checked to see whether 
the number of values chosen is sufficient. 

Now that the constants have been determined the wall 
interference may be computed. The induced velocity due to 
the wall in the direction y at right angles to the major 
axis x is 



= - P x - 3x s p 3 - 5x 4 B s - ..♦ (8) 

8=0 

Take the average between the limits x = ±s 

s 

V - ^/wdx= ^(B 1 + |k 2 B 3+ ik 4 B 5 + ...) (9) 

— s 

The value w/2 will be used in considering half of an in- 
finite horseshoe vertex filament. Furthermore, use the 
relation 

" " ' P v 2 S r 

where 

p air density 
V air velocity 




= ± p V 2 c z S 
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c 2 lift coefficient 

S area of the model 

and the upward inclination is 

La =' 



2V 



= K C i 



_1_ 
IT A 



where 

S 0 = nat, 
a 



tunnel cros s— s ect ional area 



A = (2s) /S aspect ratio 



8a 



(10) 



* = B i + T k * Bi3+ b" k * B s + • • • 



(li) 



EXACT SOLUTION FOH CLOSED AND OPEN SECTIONS 



Now, consider the more general case. By introduction 
of elliptic coordinates ( £ , n ) with the aid of the rela- 
tion z = c ch( i + in ) where c is a constant, the stream 
function for the vortex system shown in figure 2 may "be 
written as 

. T, [ch(| + |t )+eos(n-T)' )] CohU-S* )+cos(ti+ti« )] (l2 j 
Vl g [ch<| + |« )-cos(n+T)« )3 [ch(£- t» )-cos(Ti-n ' )]■- 

which may "be expanded into a convergent Eourier series for 

whore the indices 1 and ». * under the summation signs 
indicate that only odd and even terms, respectively, are 
taken* 



00 i 

^ = "tt" \ X 11-1 e~ n ^ chn| ' oo« nn !:oosnTi 
n=l 

) 



shn|' s innri 1 sinnn 



(13) 
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As in the previous case, the stream function within 
the ellipse \|/ =? + aJ/s is divided into two parts, the 
part ^2 due to the existence "of the 'boundary satisfying 
the condition A^s = 0 continuously within the field and 

at the boundary satisfying the given boundary conditions. 
The solution is generally obtained as (reference 6) 

OP 

~ n-iCAnChnlcoanti+Bjiuhn iainn.T)) (14) 

n=l 

The Fourier coefficients A n , B n are determined by the 
boundary conditions. 

In the case of the closed section the normal— velocity 
component at the boundary 



" [<] |=£o =*'^o-^'" 1/S f 

e-^^Qchn I'cosnn'sinnn 



n = l 
,52, i i 

+ 

n = 2 

00 

+ 

n=l 



\ e~~ n ^°shn s in nri 1 eosnn 

^ — A n chn| 0 sinTi + B n shn| ocosnn^ J (15) 



— i /b 

where h 2 = -^(ch 3 !- cos 2 n) must vanish. Hence, 



c 



An = - 



e 



e 



Bn shnlo 
= 0, 



(n = 


1, 


3, 


5,... 


...) . 


(n = 


2, 


4, 


6, . . . 


...) . 


(nV 


2, 


4, 


6 , . • • 


...) 


(n = 


1. 


3, 




...) 



chn| 0 



In the case of the open section the tangential velocity com- 
ponent at the boundary 
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[-^ e"" a ^°ohn|* oosnrj' cosnri 
n = l 



^,11 

OS 



e"~ n ^°shn I * sin nT| 1 sin nri 



+ ^ ^Anshn| 0 c°s an 



n = jl 

B n chn i 0 s in nr\^ J 



where h. = -ij- ( ch s | — cos 2 T|) 



-1/2 



(17) 



must vanish. Hence, 



An = 



e 



-n? 0 



s 

= 0 
e 



chnfcosnTi', (n = 1, 3, 5, ) 

(n = 2, 4, 6 , ) 



-n£, 



Bn = "chnT 

= 0 



shnl'sinnn', (n = 2, 4, 6, ) 



Yr (18) 



(n = 1, S, 5, 



) 



By determining the coefficients as in the f or egoing, the 
boundary conditions can he completely satisfied. In this 
manner the wall interference may be computed. However, 
using elliptic coordinates and obtaining as before the 
average induced velocity along the wing span is inconvenient. 
By joining the points with coordinates t, t ' and t|'' 
(n — tj 1 ) tho induced velocity in the y direction due to 
the wall is divided by the airfoil span into two parts 



OT OS 

r r— i' — I I 



TTS 



n=l 



chn £' cos nT| 



n=2 



shn I 1 s in ntj 



■] 



(16) 
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Taking the ^value w/2 for half the infinite horseshoe 
vortex and with the same notation as the foregoing. gives 

closed section: 



oo 

I 



nchnlo 



n = 1 
open section: 



ch 8 n|«cos 3 nTi 1 + ) -— ~r- sh a n| « s ln s nti ' 



n = E 



ns 



hn£ 



e - = u r 2 ! 1 ' e - n £o 

--ch^nl* cos ?nTi'~ ) ~ — r -sh a n I 1 s in 3 nT) 1 

nshnto t—> nchn? 0 

n=l n = 2 



> .(20) 



in each case: 



k £ K 

4a 



Since both of these series converge rapidly, the computation 
is simple. In the particular case that the airfoil is at the 
center of the tunnel the second series becomes unnecessary 
and the computation becomes very simple. Further, in the 
limiting case when \ ~ 00 — that is, for a circular wall — 
the foregoing relations give the results already known for 
this case. 



DEGREE OE ACCURACY OE THE APPROXIMATE METHOD 
AND COMPUTATION RESULTS 



Table I shows the high degree of accuracy obtained with 
the approximate method using the values 8 = 15°, SO 0 , 60°, 
90° in equation (7). The values in the table all refer to 
the case of the closed section with the model mounted at the 
center of the tunnel. As the ellipse approaches a circle, 
the difference is small but the law of error distribution is 
not evident* Nevertheless, an approximation with this degree 
of accuracy is sufficient for practical purposes. 

Eigures 3 to 6 all were computed from the exact equa- 
tions, the wing in all cases being mounted at the center of 
the tunnel; K gives the value of the wall interference for 
fixed aspect ratios} and 8 gives the value of the inter- 
ference for fixed ratio of wing area to tunnel area. 
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WIND TUNNEL "SECTIONS USED IN PHAOTICE 



An approximate method of solution for a closed elliptic 
wind tunnel section has been discussed in the foregoing sec- 
tion. The method, however, is of interest in the application 
to a section shape that approaches an ellipse. The tunnel 
sections used in practice consist of circular sides joined 
by their common tangents. The exact solution for this case 
would be very laborious. The approximate solution may be 
effected, however, in the following manner. As a particular 

example , take the ratio of the shorter to the longer axis 

as l*v2 and, furthermore, consider an open section. This 
corresponds to the case considered by Knight and Harris ( 
(reforenco 5), who determined the value of 8 experimentally* 

Since an open section is dealt with here instead e'f a 
closed section, the tangential instead of the normal velocity 
component at the wall must vanish. Hence, the coefficients 
B n , instead of by (?), must be determined by • 

BjCOsCS + X ) +B 3 a s cos ( 39 + X ) + B 5 a 4 cos(56 + \) + ... 

k g cos (9 + X) - a s cos(9 - X) , 
(> - k s ) 2 + 4o- 2 k £ sin s e 

where, however, a and X are not computed by the previous 
relations but for the circular part (0° <; 9 <; 67.5°) by 

a = ein(9 + X)/72"sin9, einX = ( Jz - l)sin9 

and for the straight part by 

a = 1 sfe sin6, X = 90° - 9 

Here six terms of the series are taken. By taking the value* 
9 = 0°, 15°, 30°,. 45°, 60°, 75° the coefficients in (2l) can 
be determined. Since, however, there is a discontinuity at 
6 '= 65.5°, the computation is not as simple as in the case 
of the ellipse. Figure 8 shows the curves obtained by com- 
puting the two sides of equation ( 31 ) using the value k = G,6 
and the difference between them may be seen. The difference 
as compared with the case of the ellipse is graphically in- 
distinguishable.. But the normal component at the boundary, 
as shown in the figure* increases with increasing 9, and 
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at 8 = 67.5° the ratio with the tangential component is 
about 1:23* and consequently the error in the position of 
the boundary is below 0.01a. Hence, the error in the value 
of 8 obtained by this computation may be considered to 
lie within 1 to 2 percent. Further, according to the tests 
of Knight and Harris the value of 6 obtained from the 
curve of c z against the angle of attack (denoted by 5 a ) 
differs from the value of 6 obtained from the curve of 
c z against c x (denoted by Sj) ) while according to the 

present computation the two values agree. The value of 6 
obtained from the computation is somewhat smaller than 5]), 
the computation being based on the horseshoe vortex hypotho— 
s is . 



CONCLUSION 



A problem of great interest is, for a given ratio S/S 0 
of the area of the model to the area of the wind tunnel and 
for a fixed ratio k of the span of the wing to the width 
of the tunnel, to determine the elliptical section for which 
the wall interference is a minimum. For the case of the 
closed working section near the value k = 0 the favorable 
ratio of the minor to the major diameter, as shown in fig- 
ure 3, is 1: <s/s . This agrees with the result obtained by 
Glauert for a rectangular tunnel. If k is large, however, 
this is not necessarily the case and a" flatter ellipse is 
more favorable. 

In the case of the open— section tunnel the value of 8 
for the given value of k is smallest for the circular 
tunnel. But as k becomes larger, the value of 8 for 
the elliptic tunnel gradually decreases. An interesting 
result is also the fact that for both the closed and the 
open sections the smallest 6 is obtained if the wing tips 
are at the foci of the ellipse. In the case of the tunnel 
section consisting of circular arcs joined by their cemmon 
tangents, the centers of the circles may be considered as 
the foci. 



Translation by S. Heiss, 
National Advisory Committee 
for Aeronautics, 
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TABLE 1 





J2 "b _ 1 
2 * a " 2 




e 


1 D 1 
Jz' a a/2 

V « V 




k 


K approximate 


K exact 


k . 


K approximate 


K exact 


0.2 


0.038 


0.0381 


0.2 


0.0254 


0.0254 


.4' 


.1317 


.1308 


.4 


.0969 


.0969 


.6 


.2542 


,2541 


.6 


.2110 


.2110 


.8 


.4320 


.4305 

i — . . 


.8 

- 


,3932 


.3946 



TABLE 2 



k = 


0.45 


0.60 


0.75 


5a (test) 


-0. 249 


-0.193 


-0.19 4 


6t) (test) 


-.170 


-.160 


-.164 


6 ( computat ion) 


-.148 


-.145 


-.15.0 


5 (ellipse) 


-.149 


-.147 


-.148 
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Figs. 1, 2, 3, 4, 5, 6, 7, 8 
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Figure 6.— Open eeotion 
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Figure 3.— Closed section. 
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Figure 4.— Open section. 
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